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We predict the transverse momentum distribution from proton-proton and nuclear collisions 
at the LHC. After demonstrating that higher-twist nuclear effects are very small, we propose Z° 
production as a precision test for leading-twist pQCD in the TeV energy region. We also point out 
that shadowing may result in unexpected phenomenology at the LHC. 



The physics plans for the Large Hadron Collider (LHC) 
at CERN, which is going to be the highest-energy accel- 
erator on Earth, include a heavy-ion program. Quantum 
Chromodynamics (QCD) for both hadronic and nuclear 
collision will enter a new era at the LHC, where we hope 
to discover new physics. However, "standard physics" 
needs to be tested in the new, high-energy regime. 

At LHC energies, perturbative QCD (pQCD) provides 
a powerful calculational tool |^] . Clearly, an understand- 
ing of pQCD at the hadronic collision level is a pre- 
requisite for the discussion of particle production in nu- 
clear collisions. Enhanced power corrections from mul- 
tiple scattering in both the initial and final states, not 
to mention potential new physics from the quark-gluon 
plasma (QGP), make pQCD predictions for nuclear col- 
lisions more difficult than for hadronic collisions. 

The testing of pQCD in nuclear collisions at the LHC 
requires "clean processes", where pQCD works well on 
the hadron level. One of the important recent advances 
in pQCD theory is a reorganization of perturbative cor- 
rections (folding selected logarithmic contributions at all 
orders into exact low-order calculations), which is be- 
ginning to provide practical applications The soft- 
gluon resummation for the inclusive production of color- 
less massive states [^-|| may be the best understood and 
best tested among these resummation techniques. For 
the transverse momentum distribution of heavy bosons of 
mass M, when pT <C A/, the pr distribution calculated 
order-by-order in Us in conventional fixed-order pertur- 
bation theory receives a large logarithm, ln(M^/p|,), at 
every power of as , even in the leading order in . There- 
fore, at sufficiently small p^, the convergence of the con- 
ventional perturbative expansion in powers of as is im- 
paired, and the logarithms must be resummed. 

The heavy-ion program at the LHC will make it possi- 
ble to observe the full pr spectra of heavy vector bosons 
in nuclear collisions and will provide a testing ground for 
resummation theory in nuclear collisions. In the present 
paper, we focus on production |||,|^. Based on LHC 
design luminosities Q , we estimate that a month of run- 
ning will provide ^ 4 * 10^ events in a proton-proton 
{pp) collision, and ^ 8 * 10^ Z^ events in a Pb+Pb col- 
lision in a pt interval of 0.5 GcV in the peak regions 
of the corresponding spectra. Due to the large mass of 
the and no final state rescattering in its production, 



nuclear effects from final state interactions are expected 
to be small. We will show that the power corrections 
enhanced by initial state rescattering also remain small. 
Thus, leading twist pQCD should work well here. The 
only important nuclear effect left is the nuclear modifica- 
tion of parton distribution functions (shadowing) . There- 
fore, Z'^ production could provide a bench mark test for 
pQCD at the LHC in both pp and nuclear collisions. 

Resummation of the large logarithms in QCD can be 
carried out either in p^-space directly Q, or in the so- 
called "impact parameter" , 6-space, which is a Fourier 
conjugate of the pT-space. Using the renormalization 
group equation technique, Collins and Soper improved 
the 6-space resummation to resum all logarithms as sin- 
gular as W\NP / p].) / p]. with TO > 0, when pr ^ Q 
Collins, Soper, and Sterman (CSS) derived a formalism 
for the transverse momentum distribution of vector bo- 
son production in hadronic collisions 0. In the CSS for- 
malism, non-perturbative input is needed for the large 
h region. The dependence of the pQCD results on the 
non-perturbative input is not weak if the original extrap- 
olation proposed by CSS is used. Recently, a new ex- 
trapolation scheme was proposed, based on solving the 
renormalization group equation including power correc- 
tions 1^. Using the new extrapolation formula, the de- 
pendence of the pQCD results on the non-perturbative 
input was significantly reduced. The results agree with 
Tevatron CDF jl^ and DO data very well in the en- 
tire Pt interval from pt ^ 1 GeV to pr as large as the 
the vector mass. 

For vector boson {V) production in a hadron collision 
hA + hs, the CSS resummation formalism yields Q: 

da{hA + hB^V + X) _ 1 [,27rp^i 

dA'Pdydpl (27r)2 J 

Wib,M,XA,XB)+Y{pT,M,XA,XB) , (1) 

where xa = Mj^fs and xb — e^^ M j yfs^ with rapid- 
ity y and collision energy ^/s. In Eq. (|l]), the W term 
dominates the pt distributions when pj- <C M, and the 
Y term gives corrections that are negligible for small p^, 
but become important when p^ M . 

The function W{b, M,xa,xb) can be calculated per- 
turbatively for small b, but an extrapolation to the large b 
region requiring nonperturbative input is necessary in or- 
der to complete the Fourier transform in Eq. dll). In oder 
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to improve the situation, a new form was proposed |^] 
by solving the renormahzation equation including power 
corrections. In the new formalism, W{b, M, xa, xb) = 
WP'^*(b,M,XA,XB), when b < b^ax, with 

WP^'-'ib, M, XA, xb) = e^(''*') w{b, c/b, xa,xb) , (2) 

where all large logarithms from ln(l/6^) to In(Af^) have 
been completely resummed into the exponential factor 
S(b, M), and c is a constant of order unity [Q. 

W{b, M, XA, xb) = WP'^\b,r.a.)F''P{b- braa.) , (3) 

where the nonperturbative function F^^ is given by 



F^^^exp{~HM%l^J^) gJCb')^ 
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In Eq.(||) bmax is a parameter to separate the perturba- 
tively calculated part from the non-perturbative input. 
Unlike in the original CSS formalism, W{b,AI,XA,XB) 
is not altered when b < bmax, and is independent of 
the nonperturbative parameters. In addition, the b- 
dependence in Eq. (^) is separated according to differ- 
ent physics origins. The (6^) "-dependence mimics the 
summation of the perturbatively calculable leading power 
contributions to the renormahzation group equations to 
all orders in the running coupling constant as (/i) . The b^- 
dependence of the g2 term is a direct consequence of dy- 
namical power corrections to the renormalization group 
equations and has an explicit dependence on M. The §2 
term represents the effect of the non- vanishing intrinsic 
parton transverse momentum. 

A remarkable feature of the &-space resummation 
formalism is that the resummed exponential factor 
exp[S'(6, M)] suppresses the &-integral when b is larger 
than 1/M. Therefore, it can be shown using the saddle 
point method that, for a large enough M, QCD pertur- 
bation theory is valid even at pT = ■ As discussed 
in Ref.s the value of the saddle point strongly de- 

pends on the collision energy -y/s, in addition to its well- 
known A/^ dependence. Because of the steep evolution of 
parton distributions at small x, the -y/s dependence of W 
in Eq. (|l|) significantly decreases the value of the saddle 
point and improves the predictive power of the 6-space 
resummation formalism at collider energies, in particular 
at the LHC. 

In Z'^ production, since final state interactions are 
negligible, power correction can arise only from initial 
state multiple scattering. Power corrections directly to 
the physical observables are proportional to powers of 
^QCd/Q (Q being the physical large scale). These cor- 
rections are small for production as a result of the 
large mass of the Power corrections to the evolution 
of the renormalization group equations are proportional 
to powers of AqcD/fJ-, with evolution scale /i. There- 
fore, physical observables carry the effect of the latter 



type power corrections for all /i[Qoi Q], with the bound- 
ary condition at the scale Qo- Even with large mass, Z^ 
can still carry a large effect of these power corrections. 

Equations and (Q) represent the most general form 
of VF, and thus (apart from isospin and shadowing effects, 
which will be discussed later), the only way nuclear modi- 
fications associated with scale evolution enter the W term 
is through the coefficient (72- (Since the ^2 term of Eq. (Q) 
is related to the partons' intrinsic transverse momentum, 
it should not have a strong nuclear dependence.) 

The parameters gi and a of Eq. (^ are fixed by the 
requirement of continuity of the function W{b) and its 
derivative at & = bmax- (The results are insensitive to 
changes of bmax in the interval 0.3 GeV~^ ^ bmax ^ 
0.7 GeV^i. We use bmax = 0.5 GeY-\) The value 
of 52 and g2 can be obtained by fitting the low-energy 
Drell-Yan data. These data can be fitted with about 
equal precision if the values g2 — 0.25 ± 0.05 GeV^ and 
g2 — 0.01±0.005 GeV^ are taken. As the b dependence of 
the g2 and §2 terms in Eq. (^ is identical, it is convenient 
to combine these terms and define 



G2 = ln( 1^ 
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Using the values of the parameters listed above, we get 
G2 = 0.33 ± 0.07 GeV2 for Z" production in pp colli- 
sions. The parameter G2 can be considered the only free 
parameter in the non-perturbative input in Eq. (^, aris- 
ing from the power corrections in the renormalization 
group equations. An impression about the importance of 
power corrections can be obtained by comparing results 
with the above value of G2 to those with power correc- 
tions turned off (G2 = 0). We therefore define the ratio 
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The cross sections in the above equation and in the re- 
sults presented in this paper have been integrated over ra- 
pidity (—2.4 < y < 2.4) and invariant mass squared. For 
the parton distribution functions, we use the CTEQ5M 
set in the present work. 

Figure 1 displays the differential cross sections and 
the corresponding Rq^ ratio (with the limiting values of 
G2 = 0.26 GeV2 (dashed) and G2 = 0.40 GeV^ (solid)) 
for Z° production as functions of px at ^/s = 14 TeV. 
The deviation of R02 from unity decreases rapidly as 
Pt increases, and it is smaller than one percent for both 
yfg — 5.5 TeV (not shown) and ^/s = 14 TeV in pp col- 
lisions, even when px — 0. In other words, the effect of 
power corrections is very small at the LHC for the whole 
Pt region. 
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FIG. 1. (a) Cross section da/dpT for Z° productfon in pp 
collisions at the LHC with ^/s = 14 TeV; (b) defined in 
Eq. (|) with G2 = 0.26 GeV^ (dashed) and 0.40 GeV^ (sohd). 

In lack of nuclear effects on the hard collision, the pro- 
duction of heavy vector bosons in nucleus-nucleus (AB) 
collisions should scale, compared to the production in pp 
collisions, as the number of hard collisions, AB. How- 
ever, there are several additional nuclear effects on the 
hard collision in a heavy-ion reaction. First of all, since 
the parton distribution of neutrons is different from that 
of the protons, the production cross section of heavy vec- 
tor bosons in proton-neutron interactions differs from the 
corresponding production cross section in pp collisions. 
This difference is the source of the so-called isospin ef- 
fects. At LHC, X ^ 0.02, and the magnitude of the 
isospin effects is about 2%. This is because when x is 
in this range, the u — d asymmetry is very small |l5| . 

The dynamical power corrections entering the parame- 
ter g2 should be enhanced by the nuclear size, i.e. propor- 
tional to A-^/^. Taking into account the A-dependence, 
we obtain G2 = 1.15 ± 0.35 GeV^ for Pb+Pb reactions. 
We find that with this larger value of G2, the effects of 
power corrections appear to be enhanced by a factor of 
about three from pp to Pb-fPb collisions at the LHC. 
Thus, even the enhanced power corrections remain un- 
der 1% when 3 GeV < pt < 80 GeV. This small effect is 
taken into account in the following nuclear calculations. 

Next we turn to the phenomenon of shadowing, ex- 
pected to be a function of x, the scale /z, and of the posi- 
tion in the nucleus. The latter dependence means that in 
heavy-ion collisions, shadowing should be impact param- 
eter (6) dependent. The parameterizations of shadowing 
in the literature take into account some of these effects, 
but no complete parameterization exists to date to our 
knowledge. For example, the HIJING parameterization 
includes impact parameter dependence, but does not deal 
with the scale dependence E^l^. On the other hand, 
the EKS98 |jl|l and HKM [|19[ parameterizations have 
a scale dependence, but do not consider impact parame- 



ter dependence. (The latter parameterizations have been 
compared recently |p0|.) In this paper we concentrate on 
impact-parameter integrated results, where the effect of 
the 6-dependence of shadowing is relatively unimportant 
PH , and we focus more attention on scale dependence. 
We therefore use EKS98 shadowing in this work. 
To quantify the effect of shadowing, we define 

P , , da^'''HpT,ZA/A,ZB/B) I dajpT) 

where Za and Zb are the atomic numbers and A and B 
are the mass numbers of the colliding nuclei, and the cross 
section da{pT, Za/A, Zb /B) j dpx has been averaged over 
AB, while da{pT) / dpx is the pp cross section. We have 
seen above that shadowing remains to be the only signif- 
icant effect responsible for nuclear modifications. 
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FIG. 2. Cross section da^^'^^pr, Za/A, Zb / B)/dpT for Z° 
production in Pb-|-Pb collisions at the LHC with ^/s = 5.5 
TeV, averaged over AB (solid line), compared to the pro- 
ton-proton cross section (dashed). 

Figure 2 presents da^'^'\pT, Za/ A, Zb / B) / dpr (solid 
line) compared to da{pT) / dpT in PP collisions (dashed) 
for production at = 5.5 TeV. The insert, which 
is a magnified view of the peak region of the cross sec- 
tion on a logarithmic px scale, emphasizes that the shape 
of the distribution changes from pp collisions. Most im- 
portantly, the peak moves form 3.7 GeV to 4.1 GeV. 
This small shift may be difficult to detect experimen- 
tally. However, the peak position plays an important role 
in shadowing, due to the steepness of the cross section. 
This can be seen in Fig. 3(a), which shows the shadowing 
ratio (0) (full hne) for Z° production at ^/s = 5.5 TeV. 
In Fig. 3(b) we show the Rq^ ratio defined in Eq. (||) 
for Pb-|-Pb collisions for the limiting values of G2 = 0.8 
GeV^ (dashed) and G2 = l-^ GeV^ (sohd), respectively. 
Since Fig. 3(b) provides a good measure of the overall 
uncertainty on the shadowing ratio, and this uncertainty 
is less than 2%, the characteristic shape of Rsh may be 
easier to confirm experimentally by comparing the full 
Pt spectra in Pb-|-Pb versus pp collisions at the same 
energy. 
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FIG. 3. Cross section ratios for production in Pb+Pb 
collisions at y/s = 5.5 TeV: (a) Ra^ defined in Eq. (j^ (solid 
line), and R^k with the scale fixed at 5 GeV (dashed) and 
90 GeV (dotted); (b) Rg^ defined in Eq. (|) with G2 = 0.8 
GeV^ (dashed) and 1.5 GeV^ (solid). 

The appearance of Rsh is surprising, because even at 
•pT = 90 GeV, X ^ 0.05, and we are still in the "strict 
shadowing" region. Therefore, the fact that Rsh > 1 
for 20 GeV ^ Pt ^ 70 GeV is not "anti-shadowing" ; 
rather, it is a consequence of the change of the shape 
of the cross section from pp to AB reactions shown in 
Fig. 2. To better understand the shape of the ratio 
as a function of pr, we also show R^h with the scale 
fixed at the values 5 GeV (dashed line) and 90 GeV (dot- 
ted), respectively, in Fig. 3(a). In other words, the nu- 
clear modification to the parton distribution function is 
only a function of x and flavor in the calculations repre- 
sented by the dashed and dotted lines. These two curves 
are similar in shape, but rather different from the solid 
line. In b space, W^(6, A/, a;^, xs) is almost equally sup- 
pressed in the whole h region if the fixed scale shadowing 
is used. However, with scale-dependent shadowing, the 
suppression increases as h increases, as a result of the 
scale /z ~ 1/& in the nuclear parton distribution. We can 
say that the scale dependence re-distributes the shadow- 
ing effect. In the present case, the re-distribution brings 
Rsh above unity for 20 GeV < pt < 70 GeV. When pr 
increases further, the contribution from the Y term starts 
to be important, and Rsh dips back below one to match 
the fixed order pQCD result. 

We see from Fig. 3 that the shadowing effects in the 
Pt distribution of bosons at the LHC are intimately 
related to the scale dependence of the nuclear parton dis- 
tributions, on which we have only very limited data 
Theoretical studies (such as EKS98) are based on the as- 
sumption that the nuclear parton distribution functions 
differ from the parton distributions in the free proton, 
but obey the same DGLAP evolution pSl. Therefore, 



the tranverse momentum distribution of heavy bosons at 
the LHC in Pb-|-Pb collisions can provide a further test 
of our understanding of the nuclear parton distributions. 

In summary, higher-twist nuclear effects appear to be 
negligible in Z° production at LHC energies. (The re- 
sults for production are very similar |l^].) We have 
demonstrated that the scale dependence of shadowing ef- 
fects may lead to unexpected phenomenology of shadow- 
ing at these energies. Overall, the transverse momen- 
tum distributions calculated in this paper can be used 
as a precision test for leading-twist pQCD in the TeV 
energy region for both, proton-proton and nuclear colli- 
sions. We propose that measurements of spectra be 
very high priority at the LHC. 
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